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Abstract 

Using the time-dependent-local-density-approximation (TDLDA) within a 
self-consistent linear response theory, we calculate the elementary excitation 
energies and the associated inelastic light scattering spectra of a strongly cou- 
pled two-component plasma in a double quantum well system with electron 
occupation of symmetric and antisymmetric subbands. We find, consistent 
with the results of a recent experimental Raman scattering study, that the 
intersubband spin density excitations tend to merge with the single particle 
excitations (i.e. the excitonic shift decreases monotonically) as the Fermi en- 
ergy increases beyond the symmetric-antisymmetric energy gap A$as- How- 
ever, our TDLDA calculation does not show the abrupt suppresion of the 
excitonic shift seen experimentally at a finite value of the subband occupancy 

parameter r\ = Asas/^f- 
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In an interesting recent experimental Raman scattering study, Decca et at reported 
|]J the observation of the suppression of the collective intersubband spin density excitation 
(SDE) in a coupled double quantum well (DQW) structure in the high electron density (N s ) 
limit when the two low-lying DQW subbands, the so-called symmetric (S) and antisymmetric 
(AS) levels with a single-particle energy gap < E-p (where Ep is the two dimensional 

Fermi energy) separating them, are both densely occupied. In this Rapid Communication we 
provide a detailed quantitative calculation of the elementary excitation spectra of the cou- 
pled DQW system discussing in particular the recent inelastic light scattering experimental 
observations. We find excellent agreement between our theoretical results and experimental 
data for both the charge density excitation (CDE) and the SDE spectra except for the sam- 
ples with highest N s where the subband filling parameter 77 = Asas/^f ~ 0.25. Our main 
interest in this paper, following Ref. |l|], is to investigate the strong-coupling two-component 
situation 77 < 1 when the S and AS subbands are both occupied by electrons. Around a small 
critical value of rj (^ 0.1), the experimental SDE abruptly merges with the single particle 
excitation (SPE) spectra, whereas our calculated SDE merges with the SPE monotonically 
as a continuous function of 77 as 77 approaches zero without showing any sudden collapse 
around 77 ~ 0.1. Aside from this important qualitative difference associated with the abrupt 
collapse of the experimental vertex correction (the energy difference between the SPE and 
the SDE arises from the excitonic vertex correction) around rj « 0.1, our theoretical results 
agree very well (within 0.5 meV) with the experimental measurements. 

A typical DQW structure used in our calculation is shown in Fig. 1, where the electron 
gas is taken to be confined in the translationally invariant x-y plane and the growth direc- 
tion is the z-axis. Following the experimental work of Decca et al. IJ we have used in our 
calculations several different DQW structures with varying iV s and well parameters so as 
to have many different values of Asas an d the occupancy parameter 77 (= Asas/^f)- Our 
calculation involves three steps. First, we carry out a self-consistent local density approx- 
imation (LDA) calculation @-^| of the DQW subband energy levels Ei and wavefunctions 
(f>i(z) (a typical example being shown in Fig. 1) using the best available exchange-correlation 
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potential || . Note that by definition Asas = E 2 — E\ is the energy difference between the 
lowest two subbands of the DQW. Then we calculate the irreducible (II) and the reducible 
(IT) polarizability functions of the electron layer confined in the DQW structure using the 
self-consistent LDA subband energies and wavefunctions in a linear response theory [0,f|||. 
The irreducible and the reducible response functions are formally connected by the matrix 
relation IT = lie -1 where e = 1 — VTI is the tensor representing the subband dielectric 
function with V as the direct Coulomb interaction. The irreducible response function II is 
connected to the bare LDA polarizability function Ho (which is just the two-dimensional 
Lindhard function using LDA energies and wavefunctions) through the vertex correction 
II = n (l + t/ xc n) _1 where U xc represents the LDA exchange-correlation induced vertex 
correction. Thus n , II, and IT are respectively the "bare" polarizability bubble (including, 
however, renormalized LDA quasiparticle energies), the vertex corrected polarizability, and 
the screened polarizability. Once Ho, II, and IT are calculated their respective poles (or, 
peaks) immediately give us the SPE, the SDE, and the CDE energies. The third step of 
our calculation involves directly obtaining the Raman scattering spectra for the SDE and 
the CDE which are given respectively by the spectral strengths Im U and Im II. (For the 
purpose of comparison we also calculate the SPE spectral strength given by Im n , which 
according to the simple linear response theory, should not be accessible to Raman scatter- 
ing experiments Q.) Some typical calculated Raman spectra are shown in Fig. 2 of this 
paper. The self-consistent linear response integral equations connecting the dynamical po- 
larizability functions U (z,z'), U(z,z'), U(z,z'), are solved in the subband representation. 
The procedure is standard and has been described in the literature (H|. We adapt this 
technique for the experimental DQW structures of Ref. |l] using the self-consistent LDA 
wavefunctions and energies. In the rest of this paper we discuss our numerical results com- 
paring critically with the experimental data [1]]. All our calculations assume the effective 
mass approximation for the GaAs — A^Gai-^ As — GaAs DQW structures and we use GaAs 
conduction band parameters [[J in our LDA calculations. The DQW well parameters and 
electron gas densities are taken from Ref. 0. 
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In Fig. 2 we show calculated results for the spectral weights of the SDE and SPE modes 
for several values of the wavevector transfer q for a sample with N$ = 6.35 x 10 n cm~ 2 and 
rj = 0.12. When confronted with Fig. 2 of Ref. |]J, this figure shows that the collapsed SDE 
and SPE experimental peaks show a wave-vector dependent lineshape that agrees better 
with the theoretical lineshape of the SPE than with that of the SDE. Aside from the results 
shown in Fig. 2, we only discuss results corresponding to the experimental backscattering 
geometry with very small wavevector transfer, q = 0.1 x 10 4 cm _1 as used in the experimental 
set-up |jj , in the plane of the electron layer. (All the elementary excitation energies presented 
in this paper correspond essentially to the long wavelength, i.e. q ~ 0, limit for the mode 
wavevector in the x-y plane because this is the situation reported in the work of Ref. Jl|). 
The results of this paper correspond to intersubband elementary excitations (associated with 
the quantized z-motion of the electron gas) arising from transitions between the symmetric 
and the antisymmetric energy levels (i.e. the lowest two energy levels shown in Fig. 1) — the 
elementary excitations associated with transitions to the higher subbands are at considerably 
higher energies and were not studied in Ref. |]J. Thus, the SPE peak corresponding to the 
pole in n always occurs at the energy E = Asas corresponding to the (LDA-renormalized) 
symmetric-antisymmetric gap between the lowest two levels. Note that in our calculations of 
the spectral weight functions (Im Uq, Im II, Im II corresponding respectively to SPE, SDE, 
and CDE) we use a small collisional broadening V « O.lmeV taken from the experimental 
mobility values. 

It is traditional to write the intersubband elementary excitation energies as 0,f|,|7|-|§ 

Ecb = El D + 2E SP (n s - n AS )a* (1) 



El D = E 2 SP -2E SP (n s -n AS )P* 



(2) 



where ns(AS) are the occupancies of the symmetric (antisymmetric) subbands (i.e. N$ = 
n s + n As), -^sp = A S as, and a*, (3* are parameters (which depend on E i: <t>i(z), V, and U xc ) 
which determine [10|] the depolarization shift and the vertex correction, respectively. Note 



T0|1 that the definitions of the depolarization shift (a*) and the vertex correction (/?*) shifts 
as given in Eqs. (1) and (2) above explicitly incorporate the occupancy factor dependence 
(i.e. the ns — ^as factor) arising from the Pauli principle, thus eliminating the trivial 
dependence of both the depolarization shift and the vertex correction on the occupancy 
difference as both (S and AS) subbands are occupied. Any dependence of the depolarization 
shift and the vertex correction parameters a* and (3* on the subband occupancy factor rj 
(where rj < 1 means both subbands are occupied) necessarily arises from nontrivial screening 
and exchange-correlation corrections and not as a trivial manifestation of the Pauli principle. 
Following the experimental procedure of Ref. HJ, we calculate the a* and j3* parameters by 
obtaining E$p (poles of Im U , i.e. the LDA energy levels), E^d (poles of Im II including 
vertex corrections), and Eqd (poles of Im II) from our time-dependent LDA-linear response 
calculations. (The occupancies ns and n^s are known from our LDA calculations.) 

In Figs. 3 and 4 we show a direct comparison between our theoretical calculations and 
experimental measurements for (3* and a* as functions of the filling parameter rj (Fig. 3) 
and the total electron density N$ (Fig. 4) for all seven samples employed by Decca et al. [|TJ . 
In general, (3* goes down with decrease (increase) in rj (Ns) both in experiment and theory, 
showing very good agreement (better than 0.3 meV in absolute energies of A S as — Esd) 
except for one qualitative difference, namely, in theory (3* decreases monotonically with 
decreasing (increasing) rj (N$) whereas in experiment (3* seems to go abruptly to zero around 
rj « 0.1 (i.e. the vertex correction vanishes around r\ ~ 0.1 making the SDE and SPE 
indistinguishable). Note that even around r\ ~ 0.1, the actual energetic difference between 
our theoretically predicted SDE peak and the experimental "SDE" peak (= "SPE" with 
(3* = 0) is typically small (~ 0.2-0.3 meV). For the depolarization shift (a*) experiment 
and theory, in general, agree very well (except again for rj ~ 0.25 there is some quantitative 
difference with the experimental CDE energies being typically 0.5 meV below the theoretical 
calculations). The important point to note about the depolarization shift is that a* is 
reasonably insensitive to variation in rj (or, Ns), changing little over an almost order of 
magnitude change in electron density, both in experiment and theory. 



For obvious reasons (namely, that it is not possible to vary iVs or 77 in a single DQW 
sample — each sample comes with its fixed N$ and 77 values) the experimental results (and 
the corresponding theoretical results) are for seven different DQW samples corresponding to 
the seven different data points shown in Figs. 3 and 4. Theoretically, of course, we can vary 
77 continuously in a fixed sample by changing N s continuously in our calculations. Results 
of such a calculation for a fixed sample are shown in Fig. 5 as continuous functions of Ns 
(main figures) and 77 (insets). (For completeness, in Fig. 5 we show (3* calculated using 
both local-charge-density and local-spin-density exchange-correlation potentials [[|. The 
difference between the two curves is less than the experimental errors for the range of 77 
under consideration.) These results clearly demonstrate (without the necessary numerical 
scatter of Figs. 3 and 4 which arise from using different samples) that: (i) The vertex 
correction (3* decreases monotonically with decreasing (increasing) 77 (iVg), becoming very 
small as 77 approaches zero; (ii) The calculated (3* does not agree well with experiment for 
77 ~ 0.25, and most notably, it does not go abruptly to zero for a finite value of 7/, as it does 
in the experiment for 77 pa 0.1; (iii) The depolarization shift a* is essentially insensitive to 77 
(or iVs) remaining a constant over a broad range of density and occupancy values. 

In summary, we have theoretically calculated the lowest intersubband SDE and CDE 
energies in a strongly-coupled two-component DQW structure where both the symmetric 
and the antisymmetric subbands are occupied, finding, in good agreement with a recent 
experiment 0, that the vertex correction (3* decreases monotonically with decreasing sub- 
band occupancy parameter 77 = AgAs/^-F or increasing the electron density No, whereas 
the depolarization shift a* is insensitive to changing 77 and iVg. In contrast to the exper- 
imental finding that (3* pa abruptly around 77 pa 0.1, (3* decreases monotonically in our 
theory with decreasing (increasing) 77 (N$) becoming very small for small 77. Thus, while 
our time-dependent LDA theory correctly describes the broad quantitative features of ex- 
periment quite well, the abrupt collapse of the vertex correction must be arising from higher 
order vertex diagrams (as discussed in Ref. [!|) not included in the ladder vertex corrections 
Q of our time-dependent-LDA theory. Our results for (3* can be meaningfully understood 
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as the screening out of the excitonic correction in the high density limit which makes the 
vertex correction vanish. While the excitonic correction arising from exchange interaction 
is screened, the direct Coulomb interaction leading to the depolarization shift is obviously 
unscreened, and, therefore, a* remains unchanged as N$ increases. 
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FIGURES 

FIG. 1. A typical double quantum well structure, given by the bare confining potential Vconf, 
used in the experimental study of Ref. [1], and its self-consistent LDA subband energy levels 
Ei, eigenfunctions (pi, electron density n(z), Fermi energy Ep, and effective, Hartree, and ex- 
change-correlation potentials Veff, Vh, and Vxc- The areal density is N s = 2.68 x 10 11 cmr 2 . 

FIG. 2. Calculated Raman spectra for the spin density excitation (SDE), and the single par- 
ticle excitation (SPE), given by Im II and Im n , respectively, for several values of the in-plane 
wavevector transfer q. The DQW sample has well width d w = 139A, barrier width = 28A, 
density iV~s = 6.35 x 10 n cm~ 2 and subband occupancy parameter rj = 0.12. 

FIG. 3. Calculated and experimentally measured depolarization shift and excitonic vertex cor- 
rection parameters a* and (3* as functions of the filling parameter r] = Asas/^f, for the different 
DQW samples studied in the experiments of Ref. [1]. 

FIG. 4. Calculated and experimentally measured depolarization shift and excitonic vertex cor- 
rection parameters a* and [3* as functions of the total electron density iV~S) for the different DQW 
samples studied in the experiments of Ref. [1]. 

FIG. 5. Calculated depolarization shift and excitonic vertex correction parameters a* and (5* 
as functions of areal density Ns, for a DQW structure of well width 139A and barrier width 
40A. Insets: a* and (3* for the same DQW as functions of the filling parameter r] = Asas/^f- 
The solid and dashed lines give (3* calculated using local-charge-density and local-spin-density 
exchange-correlation potentials, respectively. 
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